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Linear Collider Alignment and Survey (LiCAS) R&D group is proposing a novel auto- 
mated metrology instrument dedicated to align and monitor the mechanical stability 
of a future linear high energy e + e~ collider. LiCAS uses Laser Straightness Monitors 
(LSM) and Frequency Scanning Interferometry (FSI) [21 [3] for straightness and abso- 
lute distance measurements, respectively. This paper presents detailed simulations of 
a LiCAS system operating inside a Rapid Tunnel Reference Surveyor (RTRS train). 
With the proposed design it is feasible to achieve the required vertical accuracy of the 
order of 0(200) fim over 600 m tunnel sections meeting the specification for the TESLA 
collider g]. 

1 Principle of the LICAS-RTRS train operation 

In figure [T] the schematic view of the LiCAS train operating in the accelerator tunnel is 
presented. The train is composed of 6 cars, the distance between the centres of neighbouring 
cars is ~ 4.5 m. Each car is equipped with 4 CCD cameras and two beam splitters (BS) 
constituting the straightness monitor. The straightness monitor measures the transverse 
translation (T x , T y ) and transverse rotation (R x , R y ) with respect to a z axis defined by the 
laser beam passing through all cars in a vacuum pipe. The laser beam is reflected back 
using the retro-reflector (RR) located in the last car, illuminating the upper CCD cameras 
of the straightness monitors. 6 FSI lines placed in the same vacuum pipe between each pair 
of cars are responsible for the distance measurement along the z axis (T z ). In addition a 
clinometer located on each car provides a measurement of rotation around the z axis (R z )- 
When the train stops in front of the wall markers it firstly measures the relative position and 
rotation of all cars with respect to the first car. This defines the local reference frame of the 
train in which the location of the wall mounted reference markers are measured next. This 
procedure is repeated for each train stop. Each marker is measured up to 6 times. Finally 
the coordinates of each marker, expressed in the local train frames are transformed to the 
frame of the first train (the global frame) by fitting them to each other under the constraint 
that wall markers have not moved during the entire measurement. 



2 Opto-geometrical model of the LICAS-RTRS train 

In order to study the expected precision on the position reconstruction of the tunnel ref- 
erence markers a simulation of the LiCAS survey train was performed. To describe the 
sensing parts of the train the Simulgeo [5] package was used which allows for modelling 
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Figure 1 : Principle of the LiC AS train operation. Top view of the two train stops along the 
accelerator tunnel are presented. For detailed description see text. 

of the opto-geometrical systems. This software is also capable of performing the full er- 
ror propagation including correlations between various sub-components linked via common 
mechanical supports. 

2.1 Results of the train simulations 

Presented results were obtained assuming the intrinsic resolution of the CCD cameras and 
FSI lines equal to occd = &fsi = 1 /im. The assumed precision of the clinometer was 
&tilt = 1 fJ-rad. The simulation was performed under the assumption that all calibration 
constants (positions and rotations of CCD cameras, beam splitters, FSI light sources and 
retro-reflectors) are known to the accuracy of a pos = 1 /im for positions and a an g = 1 ^rad 
for angles. 

The long-distance operation of the train inside the accelerator tunnel was simulated by 
a set of many identical trains displaced by 4.5 m (distance between stops), each pair of them 
coupled via 5 overlapping wall markers. Simulgeo calculations provide very precise results 
(taking into account correlations between subcomponents of the system) based on the exact 
opto-geometrical model of the survey procedure. However, from the numerical point of 
view, such an approach, manipulating large matrices, is very time and memory consuming. 
The 20 train stop results (90 m tunnel section) were obtained after 34 hours of CPU time 
using 1 GB RAM memory on a 2 GHz machine (the rank of the used matrix was of the 
order of 10 000). The numerical complexity of these calculations scale like N 2 , where N is 
the number of involved coordinates. The simulation of the full 600 m tunnel section would 
require more then 7 weeks of CPU time. 
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Figure 2: LiCAS train simulation (left plot): insertion contains the results of the exact 
simulation obtained with the SiMULGEO, extrapolated on the main plot to the 600 m tunnel 
section using the formula from the random walk model. Examples of the random walk tra- 
jectories (upper middle plot) and straight line fit to selected trajectory (lower middle plot). 
RMS of the residua distribution from the straight line fit to the random walk trajectories: 
horizontal (upper right plot) and vertical direction (lower right plot). 

2.2 Random walk model 

To overcome the above mentioned limitations a simplified analytical formula inspired by a 
random walk model was derived to extrapolate the SiMULGEO predictions over long tunnel 
sections: 

f , n(n + l)(2n + l) ' ~ n{n + 1) / 9 n(n + l) 
<Jx V ,n = \[l 2 el g +o 2 xy , a z . n = ] Jcr 2 z (1) 

where n is the wall marker number, I is the effective length of the ruler (here: distance 
between cars), and the corresponding errors are the parameters of the random walk: a a 
is the angular error, a xy are the transverse errors and a z is the longitudinal error. In this 
approach the procedure of accelerator alignment resembles the construction of a long straight 
line using short ruler. The overall error is a convolution of the precision of the ruler and 
the precision of the placement of the ruler with respect to the previous measurement. The 
asymptotic behaviour of the formulae from equation no. [1] is: o xy , n ~ n~5 , and o z%n ~ n. 
This fast growth of errors (especially for transverse directions) is a consequence of the fact 
that the errors are highly correlated and the precision of the n th element depends on the 
precision of all previous points. Formulae Q] were fitted to the Simulgeo points determining 
a„, a xy , a z and then extrapolated over a 600 m tunnel section (fig.[2]left plot). The obtained 
predictions refer to the precision of the placement of the n th accelerator component with 
respect to the first one. However this is not the ultimate measure of the quality of the 
accelerator alignment. The relevant parameter is the mean deviation of each component 
from the ideal straight line which can be expected from the above procedure. To obtain the 
final prediction on the deviation of the alignment from the straight line a series of random 
walk trajectories was generated using the parameters fitted to the Simulgeo points (fig. 
middle column). A straight line was fitted to each trajectory and the corresponding residua 
were calculated. The extracted RMS values of the residua distributions for each marker 
along 600 m provide the measure of the accuracy of the whole procedure. Because of high 
correlation between errors for n th and (n + l) th marker the generated trajectories exhibit 



LCWS/ILC2007 



Mean Fourier Spectra of Random Walk Trajectories 




wave length (m) 



Figure 3: Fast Fourier Transformation (FFT) of random walk alignment trajectories. Mean 
value for position amplitude (upper plot) and residua amplitude (lower plot) from several 
Monte Carlo generations is presented. 

much smaller oscillations that would be expected from completely random process. Figure 
[2] (right column) summarises the results obtained in this analysis demonstrating that the 
vertical precision of the order of 0(100 fj,m) over 600 m is feasible. 

2.3 Spectral analysis of the alignment trajectories 

In order to study the spectra of the alignment trajectories obtained form the LiCAS Random 
Walk model the Fast Fourier Transformation (FFT) was performed. In figure [3] the mean 
values of the position and residua amplitude for several Monte Carlo generations is presented 
for 600 m tunnel section. The spectra are dominated by the long wave length components 
reaching the amplitude of about 50 fim at 600 m. In the future more realistic model should 
also include the white noise component from the short distance 'stake out' measurements. 
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